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(54) Imi 

(57) A method of reconstructing an image captured 
as a stream of image data, for example as input received 
from a linear sensor in unconstrained scanning, com- 
prises reconstructing the image in the form of a plurality 
of tiles. Each tile comprises a pixel grid of predetermined 
dimension representing a specific spatial region of the 



image. The tiles tessellate a rectilinear image space. 
Tiles can be created when required and compressed 
when no longer active, thus minimizing memory require- 
ments. Devices utilizing this method are provided. The 
method is especially appropriate for use in an uncon- 
strained hand scanner, but can also be applied to pan- 
oramic capture with a digital camera. 
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Description 

Technical Field 

5 The present invention relates to image processing methods and associated devices relating to the reconstruction 
of an image captured as image data. The present invention is particularly applicable to images captured by a scanner, 
especially a hand scanner adapted for movement across a surface in an unconstrained manner. 

Background Art 

10 

Scanners for electronically forming an image of an original are known. Typically, the captured image provided by 
a scanner is a pixel data array that is stored in memory in a digital format. A distortion-free image requires a faithful 
mapping of the original image to the pixel data array. Scanners typically use mechanical constraint during the image 
capture process in order to maximize the likelihood of faithful mapping. Types of scanners known in the art include 

is drum scanners, flatbed scanners, two-dimensional array scanners, sheet-fed scanners and hand scanners. The 
present invention is particularly applicable to hand scanners. 

Conventional hand scanners require a userto move a linear array of electrooptical sensor elements over an original. 
The movement is generally by hand manipulation, but in some embodiments is motor driven. Array-position information 
is determined using methods such as those employed in operation of a computer "mouse." As a linear sensor array is 

20 moved, the rotation of wheels, balls or rollers that are in contact with the original is sensed, and the position information 
is determined from the mechanical details ol the rotation. In general, the surface of the mechanical element in contact 
with the original has a high coefficient of friction, e.g. rubber, so as to resist slip and skid. A cylindrical roller or two 
wheels connected by a rigid axle may be used to enforce a singls translational degree of freedom during the scanning 
process. A straight-edge or other fixture is often used to fix the scan direction with respect to the original and to further 

2S enforce the translational constraint provided by the pair of wheels or the roller. Nevertheless, the position encoder 
approach is one that is often susceptible to slips and skips, so that the pixel data array lose s its correspondence with 
the image on the original. 

Hand scanners are typically connected directly to a separate computer for image data storage, processing, and 
use. The scanners provide feedback to the user, typically by means of green or red light emitting diodes, to maintain 
30 the appropriate speed for the desired image resolution. Some hand scanners use electromagnetic brakes to prevent 
the userf rom dragging the scanner over the image too rapidly, with the mechanical resistance increasing with increases 
in scanning speed. 

Hand scanners utilize relatively small imaging arrays and generally cannot handle larger than A6 documents in a 
single pass. This requires stitching algorithms to join together multiple swaths of a larger document. Swath stitching 

35 may be done in a separate operation by a separate computer. Scanning a multi-page business document or report 
with a hand scanner is a tedious process that often yields low-quality results. Techniques for stitching image swaths 
are known in the scanning art. These techniques typically require a pair of complete image swaths and produce a 
single, global transformation which brings the two swaths into registration. Improved techniques are taught in the In- 
ternational Patent Application Publication No. WO g6/27257. 

40 As previously noted, some type of fixture is typically used with a hand scanner. In the absence of a fixture, there 
is a tendency to impose some rotation as the hand scanner is moved across an original. It the user's elbow is resting 
on a flat surface during movement of the scanner, the rotation is likely to have a radius defined by the distance between 
the scanner and the user's elbow. As a consequence, the scanned electronic image will be distorted. Other curvilinear 
movements during a swath of the scanner will also create distortions. 

4S In forming a scanned electronic image it is necessary to acquire navigation information that is acquired along with 

image data. A hand scanner in which navigation information is obtained without compromise to ease of use of the 
scanner is described in US Patent No. 5578813. The navigation information is acquired by means of at least one 
navigation sensor that detects inherent structure-related properties of an original being scanned. Movement of an 
image sensor along the original is tracked by monitoring variations of the inherent structure-related properties as the 

so image sensor is moved relative to the original. The inherent structure-related properties that are monitored are inherent 
structural features, such as paper fibers, or other constituents of the original. Navigation can alternatively be speckle- 
based, wherein movement of the image sensor along the original is tracked by monitoring variations of speckle patterns 
produced using coherent illumination for acquiring the navigation information. Movement of the scanner across the 
surface of the image to be scanned can thus be essentially arbitrary, as it is not constrained. 

55 "Inherent structure-related properties" are properties of the original that are attributable to factors that are inde- 
pendent of forming image data and/or of systematic registration data on the original. The navigation information may 
be formed by generating a position signal that is responsive to detection of inherent structure-related properties, such 
as a position signal of speckle information or a position signal that permits tracking of individual inherent structural 
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features. "Inherent structural features" are those features of an original that are characteristic of processes of forming 
the original and are independent of forming image data and/or systematic registration data on the original. For example, 
if the original recorded media is a paper product, the inherent structural features of interest may be paper fibers. As 
another example, navigation of the image sensor across a glossy original or an overhead transparency film may be 

s determined by tracking surface texture variations that affect specular fields. Typically, the inherent structural features 
are microscopic, e.g. between 1 0 and 40 u,m, features of surface texture. 

The present invention is directly applicable to hand scanners as taught in International Patent Application Publi- 
cation No. WO 96/27257 and US Patent No. 557881 3, both of which applications are incorporated herein by reference 
to the extent that this is permissible by national law. 

10 Processing problems are caused by the unconstrained nature ol input in many of the embodiments set out in the 

above prior art documents. In order to be sure that the pixel array for the reconstructed image will be capable of 
representing all the data captured, it is necessary for it to be relatively large. As data capture could begin at any point 
on the image, in conventional structures it would be necessary to use a pixel grid several times larger than the image 
to be captured. It is necessary with such an approach to allocate considerable memory to storage of pixel array data, 

»s thus significantly increasing expense. 

The present invention aims to provide a method that significantly reduces the memory requirements for storing an 
image captured in such an unconstrained manner. In particular aspects, the invention provides for reduction of the 
computational cost and increases the speed of forming a reconstructed image from an arbitrarily obtained captured 
image. The method is particularly directed to freehand scanning, and freehand scanning devices. 

20 

Summary of the Invention 

According to the present invention we provide a method of reconstructing an image from a plurality of image 
elements, each image element comprising image data captured in an image sensor and position data indicating a 

25 position of that image element with respect to others of the image elements, the method comprising reconstructing the 
image in the form of a plurality of tiles, each tile comprising a pixel grid of predetermined dimension representing^ 
specific spatial region of the image. 

This position data may be derived from navigation means which provides the position data for an image element 
on capture of that image element by an image sensor. This is particularly advantageous where a linear sensor is used 

30 to capture data. In other aspects of the invention, the position data is derived from a preceding step of determining 
from features of the image elements a relative position for each image element with respect to others of the image 
elements. This is particularly advantageous where data is captured with an area sensor, and it becomes possible to 
assess relative positions of captured data by evaluation ot overlapping areas. 

Advantageously, the tiles are contiguous, do not overlap, and are adapted to tessellate in two dimensions. The 

3S method is particularly appropriate where the image is captured in the form of a stream of sensor readings, each sensor 
reading comprising image data and corresponding position data. Such data capture occurs in the case of freehand 
scanning with a linear sensor scanned across an image surface. n ' 

The method is especially appropriate to unconstrained input obtained from an arbitrary scanning path across the 
image. This is because tiles may be created at any position required, but need not be created until they are actually 

*o required. In a preferred embodiment, the method comprises the steps of: managing tiles by determining from the 
position data if tiles exist corresponding to the spatial positions of a group ol input sensor readings in Ihe image, and 
by creating one or more new tiles as necessary so that tiles corresponding to the spatial positions of the group exist; 
and writing to tiles by assigning pixel values to the pixels of tiles corresponding to the spatial positions of the input 
sensor readings of the group on the basis of the image data and the position data; wherein said step s are repeated 

45 for successive groups of sensor readings. 

In preferred embodiments, the method includes creating a tile index indexing each tile created with its specific 
spatial region of the image. Advantageously, the tile index is provided as a tree, with each node of the tree pointing to 
a plurality of locations at a multiple of the spatial resolution of that node. 

In a further advantageous aspect of the invention, tiles which have been written to can be compressed while image 

so capture is still proceeding. Advantageously, in the tile managing step the tiles corresponding to the spatial positions of 
a previous group operated on are compared with the tiles corresponding to the spatial positions of a present group, 
and wherein newly required tiles are created and wherein tiles of the previous group not p resent in the present group 
are compressed. The tiles which are compressed are written to a compressed tile storage area and are deleted from 
an active tile buffer holding tiles of the present group. 

55 The invention further comprises scanning devices adapted to carry out methods as set out above. When the present 
invention is implemented in a hand scanner which has sufficient computing power to do all of the processing, the 
method comprises using error estimates to feedback corrections to the navigation device. An exemplary form of these 
error estimates and their application is discussed in detail below with respect to an embodiment of the invention. Thus, 
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in real time, errors are identified and used to correct the navigation device itself as well as to correcl errors in the image 
being reconstructed. Alternatively, the present invention may be implemented in a scanning device for collecting image 
data and which is designed for interconnection to a computer system for performing a method according to the present 
invention. 

s A scanning device according to the invention will be described. The device may have the computing power to 

perform the processing required by a method according to the present invention or, alternatively, the device may be 
designed for interconnection to a computer system for performing a method according to the present invention. 

In the embodiment to be described, the reconstruction of a captured image includes using navigation information 
that is acquired along with image data, and then rectifying the image data based upon the navigation and image infor- 

10 mation. Rectification is a process of using navigation data to arrange acquired image data to achieve conformance 
between original and output images. In a preferred embodiment, the navigation information is acquired by means of 
at least one navigation sensor that detects inherent structure-related properties of an original being scanned. Movement 
of an image sensor along the original may be tracked by monitoring variations of the inherent structure-related prop- 
erties as the image sensor is moved relative to the original. The inherent structure-related properties that are monitored 

>s are inherent structural features, such as paper fibers, or other constituents of the original. Navigation may also be 
speckle-based, wherein movement of the image sensor along the original is tracked by monitoring variations of speckle 
patterns produced using coherent illumination for acquiring the navigation information. 

Thus, the contemplated approaches to acquiring navigation information vary in scope. In the broadest approach, 
there is no limitation to the sources of navigation information that are to be used to remove distortion artifacts of cur- 

20 vilinear and rotational movement of the scanning device along a scan path. The navigation signal may therefore be in 
the form of a position signal that is responsive to detection of image data on Ihe original (e.g., identification of edges 
of text characters), with the position signal then being used in the manipulation of an image signal. A second approach 
is one in which a posit ion signal is responsive to detection of inherent-structure related properties, such as the prop- 
erties that determine speckle patterns. The third approach is to track navigation of the scanning device by monitoring 

25 the positions of individual inherent structural features (e.g. , paper fibers) over time. 

In the embodiments to be described the image sensor is a linear array of electrooptical elements, while the navi- 
gation approach utilizes at least one two-dimensional array of navigation sensor elements. By placing a separate two- 
dimensional navigation array at each end of the image sensor, the scanner is afforded three degrees of freedom of 
movement. If the original is planar, two of the degrees of freedom are translational and are perpendicular to each other 

30 within the plane of the original, while the third degree of freedom is rotational about the normal to the plane of the 
original. The accuracy of rotation tracking is enhanced by the use ol two navigation arrays, with each array having a 
smaller array extent than would be necessary if only a single navigation array were used. While the described embod- 
iment is one in which a navigation sensor is a two-dimensional array, linear arrays may also be used. Navigation 
information for rectifying image data could feasibly be acquired by fixing other position-tracking means to the scanning 

35 device in manners independent of any inherent structure-related properties of the original. 

The navigation sensors are in a known position relative to the image sensor. Preferably, the navigation sensors 
are as close to the end points of the imaging sensor as possible, so that the navigation sensors are less susceptible 
to travelling beyond the edge of an original as the image array is moved. The image sensor forms a signal that is 
representative of an image of interest. Simultaneously, each navigation sensor forms a signal representative of the 

40 inherent structure-related properties ol the original. The scanning device may be moved in a freehand meandering 
pattern, such as one of alternating left-to-right and right-to-left movements with descent along the original, with the 
device remaining in contacl with the original. Each one of the side-to-side swaths should overlap a portion of the 
previous swath, so that the image may be manipulated with respect to position and stitched either during or following 
the scanning process. The manipulation of the image signal is a rectification of image data, with the rectification being 

« based upon the relative movement between the navigation sensor or sensors and Ihe inherent structure-related prop- 
erties detected by the navigation sensors. Stitching is a process used to connect image data acquired during successive 
swaths such that the positional relationship between data in separate swaths conforms, after stitching, to the positional 
relationship between those image data in the original image. 

Further aspects of sensor design and image processing methods with which the present invention may be em- 

50 ployed are discussed in International Patent Application Publication No. WO 96/27257 and US Patent No. 557881 3. 

Brief Description of the Drawings 

Fig. 1 is a perspective view of a hand-held scanning device following a meandering path on an original in accord- 
ss ance with the invention. 

Fig. 2 is a rearward view ol imaging and navigation sensors of the scanning device of Fig. 1 . 

Fig. 3 is a perspective view of the scanning device of Fig. 1, shown with the imaging and navigation sensors 

exposed. 
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Fig. 4 is a schematical side view of an illumination system tor one of the navigation sensors of Fig. 3. 

Fig. 5 is a side schematic view of a light emitting diode and optical elements lor providing the illumination described 

with reference to Fig. 4. 

Fig. 6 is a conceptual view of the image capture operation of the scanning device of Fig. 1. 
s Fig. 7 is an operational view of one embodiment of the navigation processing of the scanning device of Fig. 1. 

Fig. 8 is a schematic view of selected steps of Fig. 7. 
Fig. 9 is a block diagram of the components for carrying out the steps ol Fig.8. 

Fig. 10a is a representation of an increment of a position-tagged data stream typical of the output from the navi- 
gation processor of Fig. 9. 

10 Figure 10b is a representation of a buffer storing a multitude of position-tagged data increments; 

Figure 11 is a representation of a tiling scheme according to an embodiment of the invention; 

Figure 12a is a spatial representation of a tile index in accordance with the scheme of Figure 11 ; 

Figure 12b shows an index tree corresponding to the tile index of Figure 1 2a; 

Figure 13 shows an index tree as in Figure 12b represented as an array of quads; 
is Figure 14 shows a schematic representation of a method of tiling, rectifying and compressing scanned data in 

accordance with embodiments of the invention; 

Figure 1 5 shows the construction of a table of active tile indices in the context of the method illustrated in Figure 14; 
Figure 16 shows the construction of a tile cache in accordance with the method illustrated in Figure 14; 
Figure 17 shows a rectilinear image space showing the locus of the end points of the samples from the linear 
20 image sensor; 

Figure 18 shows a portion of the rectilinear image space of Figure 17 defining the scanned image with two linear 
image sensor positions mapped thereon; 

Figure 1 9 shows the backmapping of a pixel in the rectilinear space of Figure 18 on to a space defined by Ihe two 
linear image sensor positions; 

25 Figure 20 shows a portion of the rectilinear image space of Figure 17 defining the scanned image with a series of 
linear image sensor positions mapped thereon; * 
Figure 21 shows the backmapping of a pixel in the rectilinear space of Figure 20 on to a space defined by the 
linear image sensor positions; *■ 
Figure 22 shows collection of swathes by the scanning device of Figure 1 ; 

so Figure 23 shows the buffers used in a two stage process for the stitching of successive swathes; and 
Figure 24 illustrates swathes being collected when Ihe two stage process of Figure 23 is in use. 

Best Mode for Carrying Out the Invention 

3S With reference to Fig. 1 , a portable, hand-held scanning device 1 0 is shown as having followed a meandering palh 
12 along an original 14. The original may be a piece of paper, an overhead transparency film, or any other image- 
bearing surface. The inherent structure-related properties of the original may be used to provide position information 
during navigation along the meandering path. In this approach, the positions of inherent structural features are tracked 
and the position information is used to rectify image data, but other embodiments will be described. The scanning 

40 device is preferably sell-contained and battery powered, but may include a connection to an external source of power 
or to data ports of computers or networks. 

The scanning device 10 of Fig. 1 includes an image display 16 to allow viewing of a captured image. However, a 
display is not essential to the use ol the scanning device. 

Referring now to Figs. 1-3, the forward side 18 of the scanning device 10 includes a pivoting member 20 that aids 

45 in maintaining proper contact between the original 14 and an imaging sensor 22. The imaging sensor 22 comprises a 
linear array of image sensor elements. Navigation sensors 24 and 26 are located at the opposed ends of the imaging 
sensor. The navigation sensors 24,26 are mounted on the pivoting member and so the navigation sensors are in a 
fixed location relative to the imaging sensor. 

The scanning device 10 allows three degrees of freedom, with two being in translation and one in rotation. The 

so first degree is the side-to-side movement (X axis movement) along the original 14. The second degree of freedom is 
movement upwardly and downwardly along the original (Y axis movement). The third degree of freedom is the ability 
to operate the device with rotational misalignment the imaging sensor 22 relative to an edge of the original 1 4 That is, 
the imaging sensor 22 may have an angle of attack that is not perpendicular to the direction of device translation. 
For reasons of physical compactness, the imaging sensor 22 is preferably a contact image device, but for appli- 

55 cations in which compactness is less of a concern or a smaller image is desired, sensors employing projection optics 
may be employed, with magnification less than unity. In such applications, the elements of the imaging sensor 22 
should be smaller and more closely packed together. Contact imaging devices typically employ lenses sold under the 
trademark SELFOC, which is a trade mark of Nippon Sheet Glass Company Limited. Less conventionally, contact 
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imaging can be obtained using interleaved array elements of sources and proximal sensors, without any imaging lenses. 
Conventional imaging sensors for scanning applications may be used The imaging sensor may be part of a unit that 
also includes an illumination source, illumination optics, and image transfer optics. 

The imaging sensor 22 is shown as a linear array of discrete optically sensitive elements. The spacing of the 
elements plays a role in determining the spatial resolution of the scanner 10. For example, a linear array having a 
length of 101 .6 mm requires 1200 sensor elements to achieve a resolution of 300 dpi. The sensor may be a charged 
coupled device, an amorphous silicon photodiode array, or any other type of linear array sensor known in the art. 

A key consideration in the design of the imaging sensor unit is speed. The imaging sensor 22 preferably is able 
to image each pixel at approximately 5K samples per second. Linear imaging arrays generally produce a serial data 
stream, wherein pixel values, i.e. charges, are placed into a shift register and then shifted out. Achieving the desired 
speed requires either very fast serial transfer rates out of the entire image array or multiple taps, so that pixel values 
can be shifted through fewer cells. This introduces parallelism, which is advantageous for digital processing. 

Another consequence of the speed requirement is that the products of pixel areas, at the surface of the original, 
and their solid angles of emitted light collected and conveyed to each array element should be sufficiently large to 
produce a detectable signal in integration times on the order of 200 microseconds. An enhancement option is to add 
an optical element to the sensor to increase the effective fraction of sensor pitch for which each sensing element is 
responsive. Since there is typically unused area in the array matrix, such light collection optics increase sensitivity. 

A straightforward modification of the imaging sensor 22 enables sensing of color images. Three linear arrays 
parallel to one another, each with at least one embedded filter element that selectively passes red, green and blue 
components of the incident light respectively, would allow color imaging. Alternatively, a single array having broad- 
band sensitivity may be sequentially illuminated by red, green and blue light sources. 

Regarding illumination for improving the operation of the imaging sensor 22, a linear array of high intensity light 
emitting diodes at the amber wavelength may be used. However, the selection of the preferred illumination source and 
. any optical elements is dependent upon the medium of the original. The wave-length of the light is selected to maximize 
the contrast image data acquired during the scan of a given area of the original 14, while disregarding unwanted signals. 
Illumination optics may consist of LED dome lenses or may include a light pipe consisting of a precision-molded oplical 
element that channels the illumination onto the original with a minimal amount of light loss. Such a design can aflord 
a relatively uniform illumination of a target area of the original at a wide range of angles, but blocks normal incident 
rays in order to avoid specular surface reflections. 

In Fig. 1, the meandering path 12 is shown as having four and a fraction swaths, i.e. side-to-side passes across 
the original 14. A useful imaging sensor 22 for most opportunistic applications has a length within the range of 25.4 
mm and 101 .6 mm. If the sensor 22 has a length of 63.5 mm, an A4 paper can be scanned in four or five swaths. As 
will be explained more fully below, the swaths should include regions of overlap, so that a stitching process can be 
used to produce a faithful reproduction of the original image. 

The scanning device 10 typically includes at least one navigation sensor 24 or 26. In preferred embodiments, the 
device includes a pair of navigation sensors, with the sensors being at opposite ends of the imaging sensor 22. While 
a one-dimensional array of optoelectronic elements may be used, in the preferred embodiment, each navigation sensor 
is a two-dimensional array of elements. The navigation sensors 24 and 26 are used to track movement of the scanning 
device 10 relative to the original 14. 

In preferred embodiments, each navigation sensor 24 and 26 captures images related to inherent structure-related 
properties of the original in order to produce information related to the position of the scanning device 10. For most 
prior art scanning devices, inherent structural features are considered to be noise. For the scanning device 10 of Figs. 
1 -3, such features are noise with regard to the imaging sensor 22, but may be used to provide a ba sis for the navigation 
sensors 24 and 26 to generate position information. Useful, high-contras t images of surface texture can be generated 
by detecting structural variations that are inherent to the medium or are formed on the medium, e.g., text. For example, 
images may be formed based upon the contrast between shadows in valleys and bright spots at Ihe peaks of inherent 
structural features. Such features are typically microscopic in nature, often ranging between 1 0 urn and 40 u.m in size 
on common printed media. As an alternative, speckle may be used, since specular reflection of a coherent beam 
produces a contrast pattern of bright and dark areas. A third source of contrast information is color. Color contrast is 
independent of surface texture. Even when illuminating the texture-free surface with light in the visible range, color 
contrast exists between regions ol different colors, e.g., between different shades of gray. 

However, it is contemplated to use the invention for applications in which navigation information is independent of 
inherent structure-related properties of the original. Non-imaging approaches can also be used to acquire and process 
X, Y and angular position information. Such alternatives are discussed in US Patent No. 557813. 

The navigation sensors 24 and 26 of Fig. 2 effectively observe moving windows on the image of the original 14 
and produce an indication of the dispfccement in two planar dimensions between successive observations. As will be 
explained more fully below, pwel values from the navigation sensors are processed to determine proper mapping of 
image data from the imaging sensor 22. Values of a particular pixel and its nearest neighbors are used to produce an 
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array of correlation values at each pixel location. The correlation value s are based upon comparisons between a 
current image of the surface structure and a stored image representing a known position of inherent structural features, 
wherein the stored image serves as a position reference. However, operations other than the correla- tion process may 
be empbyed in manipulating the input image data to form the output image. 

s Referring now to Figs. 4 and 5, navigation sensor 24 is shown as being operativefy associated with illumination 

optics. If the original 14 is a paper product for which paper fibers are to be detected by the navigation sensor 24, the 
introduction of light at a grazing angle of incidence is preferred. While not essential, one or more light emitting diodes 
(LEDs) 28 may be used. The grazing angle 30, which is the complement of the angle of incidence, is preferably in the 
range of zero degrees and fifteen degrees, but Ihis may change depending upon the properties of the original 14. In 

10 Fig. 5, the source 28 is shown with illumination optics 34. The optics may comprise a single element or a combination 
of lenses, filters and/or holographic elements to accomplish suitable collimated and generally uniform illumination of 
the target surface. The wavelength of the light emitted by the source 28 should be selected to enhance the spatial 
frequency information available for navigation. Fixed pattern noise in the illumination field should be minimized. The 
output of source 28 may require adjustment to accommodate wide dynamic ranges of reflectivity of the medium as the 

'5 scanning device proceeds over printed materials with absorbing or reflecting inks or other marking agents. 

In Fig. 4, light from a source 35 is collimated at illumination optics 36 and then redirected by an amplitude splitting 
beam-splitter 37. That portion of the light energy from the LED directly to and transmitted through the beam-splitter is 
not shown in Fig. 4. The light energy from the beam-splitter illuminates the original 1 4 along the normal to the surface. 
Also represented in Fig. 4 is the portion of the light energy that is reflected or scattered from the original 1 4 and 

20 passed through the beam-splitter 37 for aperturing and filtering at element 38 and focusing to an image at element 39. 
The portion of light energy passing from the original to the beam-splitter and reflecting from the beam-splitter is not 
shown. The magnification of navigation imaging optics should be constant over the field-of-view of the two-dimensional 
sensor array 24 which detects the focused light. In many applications, the modulation transfer functions, i.e. the am- 
plitude measure of optical frequency response, of the navigation optics must be such as to provide attenuation before 

25 the Nyquist frequency that is determined by the pitch of Hie sensor elements of the navigation sensor and by the 
magnification of the optical elements. The optical elements should also be designed to prevent background illumination 
from creating noise. Note that a wavefront splitting beam-splitter could also be used. 

The selection of the angle of incidence depends upon the material properties of the original. Grazing angles of 
illumination generate longer shadows and more apparent contrast, or AC signal if the surface of the original is not 

30 glossy. The DC signal level, however, increases as the illumination angle approaches the normal to the original. 

Illuminating the target region of the original 1 4 at a grazing angle 30 works well for applications in which the surface 
of the original has a high degree of unevenness at the microscopic level. For example, the introduction of light from 
the source 28 at a grazing angle provides a high signal-to-noise ratio of data related to inherent structural features 
when the original is stationery, cardboard, fabric, or human skin. On the other hand, the use of incoherent light at a 

35 normal angle of incidence may be preferred in applications in wh ich position data is needed to track scanner movement 
along such originals as photographs, glossy magazine pages, and overhead transparency films. With normal illumi- 
nation, using incoherent light, viewing the original in the specularly reflected field will provide an image that is sufficiently 
rich in texture content to allow image and correlation-based navigation. The surface of the original has a microscopic 
relief such that the surface reflects light as if the surface were a mosaic of tiles, or facets. Many of the "tiles" of an 

40 original reflect light in directions slightly perturbed from the normal. A field of view that includes the scattered light and 
the specularly reflected light can thus be modeled as though the surface were composed of many such tiles, each tilted 
somewhat differently with respect to the normal. This modeling is similar to that of WW. Barkas in an article entitled 
"Analysis of Light Scattered from a Surface of Low Gloss into Its Specular and Diffuse Components," in Proc. Phys. 
Soc, Vol. 51, pages 274-292 (1939). 

& Fig. 4 shows illumination by a source 35 of incoherent light, which is directed along the normal of the surface of 
the original 14. Fig. 5 describes illumination at a grazing angle 30. In a third embodiment, no illumination is provided. 
Instead, the naviga- tion information is accumulated using background light, i.e. light from the environment. 

In afourth embodiment, coherent illumination is introduced at normal incidencetopermit speckle-based navigation. 
Relative motion between a scanning device and an original may be tracked by monitoring motion of speckle relative 

so to the navigation sensors. If coherent illumination is used without using imaging optics, then by selecting a small area 
of illumination and by having a relatively large separation between the surface of the original and the photodetector 
array of the navigation sensor 24, the resulting predominant speckle cell sizes with coherent illumination are sufficiently 
large to satisfy the Nyquist sampling criterion. The use of a beam splitter allows the direction of both the incident 
illumination and the detected scatter to be near to normal to the surface of the original, as similarly accomplished in 

55 Fig. 4. 

Referring now to Fig. 6, the scanner 10 is shown as being moved across an original 44 having a block 46 imprinted 
onto a surface of the original. Because the scanner 1 0 is not subjected to any kinematic constraints in the plane of the 
original, there is a tendency for a user to follow a curvilinear path across the original, as when the hand and forearm 
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the 



of ths user relate about the elbow. In Fig. 6, the scanning device is shown as following a curved path 48 across the 
tllwr 

is moved to the second position 52, shown in dashed lines. However 
The captured image 50 would be the stored image in the absence ot processing to be described I betew However, 
as the maTng sensor captures data related to the block 46, navigation information is acquire . n the tk«M» 
bo m nt o eormcrenavigationsensorscapturedatarelatedtoinherentstructuralfeaturesolt^ 

me« 

of the imaging sensor relative to the block 46. A faithful captured .mage 54 may then be formed. The image 54 

successive rames of navigation information, such as data related to inherent structural features The cor elations 
positions of I inherent structural features in successive frames to prov.e £™J»^ b £ 
position of a navigation sensor at a particular time. The navigation mformation ,s then used to rectify image data, 
orocessino of Fiq 7 is typically performed for each navigation sensor. 

In a 56 a reference frame is acquired In effect, the reference frame js a start posrtrcn. The posrtion of a 

navX« 

sensor atthelatertimeandthen computing correlations 60 between the reference frame and the later-acquired sample 

'^Acquiring the initial reference frame 56 may take place upon initiation of the imagmg process For example, the 
acqiLmay be triggered by mere P .acement o, the scanning device into ^^^^^^ 
scanning device may include a start button that initiates the ,mage process and the T^^J^JJS'S 
also take olace bv a periodic pulsing of the illumination system of each navigator. If there is a reflected signal that 
Zl ^l J^ of reflection or a correlation signal that indicates motion, the reference frame ,s then 

Mhough the navigation processing is performed computa.iona.ly, .he concepts of this embodim ent may be de- 
scribed wrth reference to Figs. 7 and 8. A reference frame 62 b shown as having an ,mage of a T-shaped nheren 
«aCure64. Thesizeof the reference frame ^^^^^'^J^^Z 
hescanninqdevice thedominantspatialfrequenciesintheimagingofthestructuralfeatures, and the W resolution 
^zo thi reeto^m. f «m. for « nMigMlan mmot M C thirty*^ p»l. (N) by «d^1our 

PiXe At?!aSS 

but wiic "s b tantially the same inherent structural feature, The duration dt is preferably 
IZlp^nX of the T-shaped feature 64 is less than one P «e. of the nav.ga ion ^ 
translation of Rescanning device. An acceptable time period is 50 us for velocrties of 0.45 meters/sec at 600 dp.. Th.s 
relative displacement is referred to herein as a "microstep." rD<oronrP , rame S 2 and acouirinq 

Ifthescanningdevicehasmovedduringthetimepenodb^ 
58 the sample frame 66, the first and second images of the T-shaped feature will be ones in which the feature has 
i« 

schematic representation of Fig. 8 is one in which the feature 64 is allowed to shift up and to .he nght by one p«el. 
ThP fiiii-nixel shift is assumed here only to simplify the representation. 

FJemeTt 70 in Fig 8 represents a sequential shifting of Ihe pixel values of frame 68 into the eight nearest-neighbor 

pixelflais^ 

s f etc m hi manner, the pixel-shifted frames can be combined with the sample frame 66 to produce te array 72 
o 1 pos o f ameT The position frame designated as "Position 0" does not include a shift, so that the result is merel 
IrZTl Z of frames 66 and 68 "Position 3" has the minimum number of shaded pixels, and therefore ,s the Irame 
w « 

acquired reference frame 62, which implies that the scanning device has moved leftwardly and downwardly during 

^Me other correlation approaches may be employed, an acceptable approach is a 'sum of the squared d^e, 
ences' correlation For the embodiment of Fig. 8, there are nine correlation coefficients (C k - C 0 , C, ... L a ) tormeo 
from the nine offsets at element 70, with the correlation coefficients being determined by equation. 
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where Sy denotes the navigation sensor-measured value at the position ij of the sample frame 66 and R| ( denotes the 
navigation sensor-measured value at the frame 68 as shifted at the element 70 in the k direction, with k being the 
identifier of the shift at element 70. In Fig. 8, k=3 provides the correlation coefficient with the lowest value. 

Correlations are used to find the locatbns of identical features in successive frames in order to determine the 

s displacements of the features from frame-to-frame. Summing or integrating these displacements and correcting for 
scale factors introduced through the design of the relevant optics determine the displacements of the imaging sensor 
as a scanning procedure progresses. 

As previously noted, the frame-to-frame correlations are referred to as 'microsteps,' since frame rates are chosen 
to be sufficiently high to ensure that the displacements do not exceed the dimension of a single pixel. Oversampling 

io can provide sub-pixel displacement precision. Referring to Fig. 7, a determination 74 ol whether a microstep is to be 
taken is made following each compulation 64 of the correlations. If a microstep is required, the reference frame is 
shifted at 76. In this step, the sample frame 66 of Fig. 8 becomes the reference frame and a new sample frame is 
acquired. The correlation computation is then repeated. 

While the process provides a high degree of correlation match, any errors that do occur will accumulate with each 

'5 successive shift 76 of a sample frame 66 to the reference Irame designation. In order to place a restriction on the 
growth rate of this "random walk" error, a sample frame is stored in a separate buffer memory. This separately stored 
sample frame becomes a new reference frame for a subsequent series of correlation computations. The latter corre- 
lation is referred to as a "macrostep." 

By using macrosteps, a more precise determination of scanner displacement across a distance of m image frame 

20 displacements, i.e. m microsteps, can be obtained. The error in one macrostep is a result of a single correlation cal- 
culation, whereas the equivalent error of m microsteps ism 14 times the error in a single microstep. Although the average 
of errors in m microsteps approaches zero as m increases, the standard deviation in the average of errors grows as 
m* 5 Thus, it is advantageous to reduce the standard deviation of accumulated error by using macrosteps having m as 
large as practical, as long as the two frames that define a macrostep are not so far spaced from one another that they 

2S have no significant region of common image content 

The sampling period dt does not have to be constant. The sampling period may be determined as a function of 
previous measurements. One method that employs a variable dt is to improve the accuracy of displacement calculation 
by keeping the relative displacement between successive reference frames within certain bounds. For example, the 
upper bound may be one-pixel displacement, while the lower bound is determined by numerical roundoff considerations 

30 in the processing of the navigation data. 

Referring to Fig. 9a, the image signal generated at the imaging sensor 22 may then be 'position tagged" based 
upon the navigation data. In one embodiment, pixel values from the navigation sensors 24 and 26 are received by a 
navigation processor 80 for performing the operations of Figs. 7 and 8. Based upon the computed correlations, co- 
ordinates are determined for the current position of the first navigation sensor 24 and the second navigation sensor 26: 

35 Within a macrostep the navigation processor 80 directly recovers the translational component of the motion of 
each navigation sensor. The data from both navig ation sensors must be integrated to obtain an absolute position 
estimate that also takes into account any rotation of the scan head. Processor 80 maintains a model of the current 
orientation of the scan head with respect to its initial orientation. The individual eslimates of the translation of each 
navigation sensor are interpreted with respect to this model. In turn the orientation of the scanner model is itself peri- 

40 odically updated. 

In one embodiment the macrosteps of each navigation sensor are synchronised so that if the navigation sensor 
processor 80 requires a macro step in one navigation sensor, it also triggers a macrostep in the other. This simplifies 
the interpretation of the translations recorded by each navigation sensor. If at macrostep T the orientation of the scanner 
is 0 degrees from vertical, then the recorded translation within the microstep of each navigation sensor is interpreted 
45 as follows. 

Figure 9b shows the motion of the navigation sensors as a pair. For each we have recorded translations with 
respect to a co-ordinate frame (with unit vectors labelled u and v) that describes orientation of the scanner. These have 
magnitude (u1 , v1 ) and (u2, v2) for the first and second navigation sensors respectively. The objective is to interpret 
these translations to give the updated locations of the navigation sensors with respect to the global co-ordinate frame, 
so in essence we are reconstructing the path from a piece-wise linear approximation of it, where the magnitude of the 
linear sections is determined by the size of the macrostep. 

Each unit vector is related to the orientation of the scanner 0 according to the following standard trigonometric 
relationships 
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u = (cose, -sine) 

but as we shall see estimates of u and v can be maintained without computationally expensive trigonometric calculation. 
5 In Figure 9b, P1 and P2 represent the navigation sensor locations, with respect to the global frame, at the start of 

the macrostep; they have values (x1, y1) and (x2, y2) respectively. The updated locations, with respect to the same 
global frame, some navigation sensor frames later are given by P1 ' and P2', which have values (x1 ', y 1 ') and (x2\ y2') 
respectively. 

In the limit (which we assume we are close toin orderto perform accurate integration) the v component of translation 
10 of each navigation sensor must be the same (shear affects each end equally, and in the limit rotation induces no v 
component at all). Hence we can set v1 and v2 to their mean value 



v=(v1 +v2)/2 

In which case the updated endpoint locations are given by 



P1' = P1 + u1u + vvandP2' = P2 + u2u + vv 



It is also 
of the image 



. convenient at this point to translate the position estimates for correspondence to the physical endpomts 
of the image sensor rather than the individual navigation sensors. This is done using calibration data that relates the 
physical locations of the navigation and image sensors. For simplicity, the end point positions are made to correspond 
to the centres of the first and last sensor elements of the image sensor. 

Periodically we must update the inertial frame. This should only be done at the end of one macro frame or at the 
start of the next. It requires that the macrosteps be synchronised: that is if one navigation sensor triggers a new mac- 
rostep, because it has exceeded the maximum number of microsteps, a macrostep should also be triggered in the 
other navigation sensor even if has remained stationary throughout. This approach has the additional virtue that the 
rotation of the navigation frame is always small as the maximum rotation experienced must be always be less than the 
maximum number of microsteps divided by the baseline between navigation sensors. 

Updating the inertial frame has the effect of adding back the lost rotational component of the motion of the scanner. 
This is determined implicitly from the differential translations experienced by the individual navigation sensors. Rather 
than explicitly calculating the rotational component of the motion, the inertial cc-ordinate frame is updated directly from 
its relationship to the principal axis of the scanner (the line joining the navigation sensor locations). The v axis of the 
inertial frame is a unit vector along the principal axis going from navigation sensor 1 to navigation sensor 2 and the u 
axis is the unit vector orthogonal to it that completes the 2D basis. 
The updated value for v is given by 

v= (V' , v j = (P2 , -pr , )/iP2'-pri = (pe-pryo 



and the orthogonal vector u by 



tf = (v' y , -v' x ) 

The length IP2-P1 1 is the constant (modulo navigation errors) length of the scanner, labelled D in Figure 9b. This 
greatly simplifies the calculation of the inertial frame as it avoids need for calculating the square root of the sum of 
squares for each macrostep. 

The navigation processor 80 also receives N pixel values of the imaging sensor 22 via a pixel amplifier B2 and an 
analog-to-digital converter 84. Although Fig. 9 shows only a single tap from the image sensor 22 and a single NO 
converter 84, multiple taps, each with an A/D converter, are within the scope of the invention. The current position 
coordinates are "tagged" at the ends of a line of data that corresponds to the number of pixels within the imaging 
sensor The output 86 of the navigation processor 80 is therefore a position -tagged data stream. 

The navigation processor 80 need not necessarily receive pixel data from the imaging sensor 22 and position 
information from the navigation sensors 24,26 at the same frame rate. The stream of position data may be used to 
update the estimate of the location of the endpoints of the linear sensor, which is generally at a higher spatial resolution 
than the pixel pitch of the linear sensor array. The optimum clocking frequency of the linear image sens or will generally 
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be determined by the rate of motion of the scanner, in order that the fastest moving end is sampled at an amount just 
less than the pixel pitch of the linear sensor itself - otherwise either over-sampling or under-sampling will result. Over- 
sampling will result in increased computational and data processing requirements further along the image processing 
pipeline with little or no improvement in image quality, while under-sampling will result in a definite reduction in image 
s quality. 

One solution is to sample the linear sensor at a fixed frequency that would always result in over-sampling, that is 
it would exceed the sampling requirements of the fastest possible scanning speed. Then if the scanner is moving more 
slowly than the maximum scanning speed, unwanted samples are discarded by the navigation processor on the basis 
of the current endpoint position estimates. That is, the magnitude of the relative position of each endpoint since the 

10 last time a valid increment of linear sensor samples and tagged position data was output from the navigation processor 
80 is used to gate the output of further increments. The simplest approach is to only allow output when one or other 
magnitude strictly equals or exceeds the pixel pitch (or a significant fraction of it). This approach is likely to result in 
under-sampling, particularly when scanning at just less than maximum scan speed. Alternatively, in order to prevent 
under-sampling, the current sample could be output if a prediction of the next linear sensor sample's relative end point 

'£ positions exceeds the pixel pitch, where the prediction is based on the rate of change (or higher order numerical 
derivatives) of end-point position. A third approach which also prevents under-sampling is to introduce an increment 
buffer into the navigation processor, 80, so that the previous increment can be transmitted if the relative position data 
of the current increment exceeds the pixel pitch. 

Each of the above sampling strategies are based upon synchronous sampling and will in general result in some 

20 degree of over or under sampling. A better overall solution that more closely achieves an ideal spatial sampling is for 
the navigation processor 80 to trigger the linear sensor capture asynch ronously. One approach is for the relative end 
point position since the last valid capture (or a future prediction thereof) to be used directly to identify the best point in 
time to trigger the capture and sampling of the linear sensor. 

In Fig. 10a an increment 88 ol the data stream is shown as having position coordinate cells 90, 92, 94 and 9S at 

2S the opposite ends of N pixel cells, although this ordering is not essential. 

The imaging sensor 22 is clocked as the scanning device moves across an original. As is discussed above, the* 
clocking ensures that the fastest moving element of the sensor samples at least once per pixel displacement. 

The position-tagged data stream at the output 86 of the navigation processor 80 is buffered as shown in Fig 10b. 
A buffer B holds a multitude of the increments 88 of Fig 10a. The buffer B may be used to accommodate an entire 

30 scan in which case data compression may be used. This data is than mapped to rectilinear co-ordinates to form a final 
rectified reconstructed image in accordance with a method according to the invention, as will now be described. 

In the prior art methods taught in International Patent Application Publication No. WO 96/27257 and US Patent 
No. 5578B13, the mapping to rectilinear co-ordinates is achieved by mapping to a rectilinear grid of predetermined 
size. Representing the rectilinear image as a single 2 dimensional array in this way has the twin disadvantages that it 

35 requires a large amount of working memory and restricts the scan path to fit within a pre-defined page. Accordingly, 
embodiments of the invention provide an improved method employing a dynamic tiling system. In such a system, the 
image space is represented as a number of discrete tiles which are allocated on a need to have basis at run time. '" 
With such atile based representation the actual path of data capture (a scan, in the case ol the scanner embodiment 
described above) does not need to be limited ahead of time. Using a generic tile indexing scheme as described in 

40 detail below, it is possible with a hand scanner such as described above to scan any shaped region: the overall dimen- 
sions of the scan are limited only by the total area under the scan exceeding the available memory. Furthermore, in 
embodiments of the invention, when tiles become "inactive" and are no bnger required to determine the processing 
of other tiles (i.e. they are no longer required for the purpose of rectification), these tiles can be compressed "on the 
fly" while data capture is continuing. This allows the total working memory required to store the rectilinear image during 

45 processing to be greatly reduced. 

In a preferred method each active tile is represented as 2 dimensional array of bytes, where each byte represents 
a pixel. Other pixel representations, with either multiple pixels represented as a single byte or multiple bytes per pixel, 
could easily be accommodated. For convenience, tiles are chosen to be square with a dimension (height and width) 
that is a power of 2. 64x64 is a particularly appropriate choice: the embodiments discussed below employ this dimen- 

so sion. The rectilinear image location associated with the top left hand corner of each tile is from the geometry of the 
scheme a multiple of the tile size. Tiles do not overlap; they are contiguous with respect to the rectilinear image they 
represent. 

Individual tiles exist in one of three states: 

55 empty implying that the tile does not currently exist 

active a currently active tile that is represented by a 2 dimensional array 
compressed a compleled tile that has already been compressed 
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Indexing of tiles is important for effective operation of such a method. In this regard, the approach taken by a 
preferred tiling scheme has two main components: 

♦ an indexing scheme that maintains the list of tiles that represents the scanned image and provides a method to 
5 index individual tiles spatially; and 

• a fast local index into currently active tiles that minimises the overhead of using the tile based representation for 
the frequent access into the image that occur during the rectification process. 

Each rectilinear image co-ordinate is represented as two unsigned 16 bits values; one for x and one for y. For the 
w case of 64x64 tiles, the top ten bits of x and y provide the tile index, which is used to identify the tile within which the 
rectilinear pixel is located. The bottom 6 bits of x and y give the within tile pixel location, and are used as indices for 
the 2 dimensional array that represents an active tile. If the tile identified by the top ten bits is empty, or has already 
been compressed, then it is not possible to use the bottom 6 bits directly to index the pixel in question. In these two 
cases memory must be allocated for the tile or the tile must be decompressed, respectively. An alternative to decom- 
»s pression is the allocation of an additional tile which must later be merged with the earlier-formed tile. Decompression 
and merging are computationally expensive operations and for efficient operation should be employed as little as 
possible. 

10 bits allows indexing of a tile space that has 1024 tiles on a side. Each 64x64 tile is roughly 1/5* of an inch on 
each side at 300dpi which produces an addressable space over 200 inches wide and tall. This virtual image size will 

so satisfy the most demanding scanning task, for example, if the origin of a scan path is located close to the centre of the 
address space, it is possible to scan in any one direction up to 1 00 inches before running out of tile space. The problems 
of exhausting working memory or long term storage are thus more significant: these are discussed in more detail below. 

A particularly convenient way to represent the tile index is as a 2 dimensional array. This form of tile index will 
require less memory than would a totally flat array to index equally si2ed flat images (ignoring the image storage itself). 

25 However, if all ten bits were used even this approach would require 2MB of working memory just for the index (allowing 
16 bits for each tile entry). For applications such as handheld scanning, use of such large amounts of memory is 
currently very disadvantageous. It is therefore desirable to find an alternative approach which consumes less working 
memory. 

It has been found that employing a smaller array that uses only 8 bits to index the tiles is possible. This requires 

30 only 1 30 kB of working memory and allows scans with a width and height of up to approximately 25 inches. The following 
description indicates how an indexing scheme can be constructed in this way. Figure 11 illustrates such an array 201 
(but with only 4 index bits shown). Each 16 bit entry in the tile index array is either set 1o zero (which implies that it is 
empty, as for array entries 202) or contains an index (those array entries 203 shown hatched in Figure 11 ) into a tile 
directory 21 1 . Using 1 6 bits as an index into the tile directory 2 1 1 allows access to up to 64K tiles, which limits the total 

35 area of the scan to at mosi 2,560 square inches for hand scanner embodiments at 300 dpi. Given that a letter size 
page is of the order of 100 square inches (2500 tiles) this does not pose a problem, so it is even possible for a smaller 
tile directory to be used in practice. 

Each Tile Directory entry contains a flag to indicate whether it has been compressed or not, a pointer to its start 
address and length parameter which gives the number of bytes used to store the compressed version. It is convenient 

40 to keep active and compressed data blocks in a separate chunks of memory. Eventually, upon completion of the scan, 
all blocks will be compressed so keeping the compressed data distinct from the working memory minimises extra 
processing required to coalesce the compressed memory and transfer it to long term storage. A pointer nextAvaila- 
bleDirectoryEntry indicates which directory entries have been used. They are allocated sequentially as each new tile 
is encountered and are each used to represent a single tile per scan. 

45 A particularly preferred approach is to represent the index as a partial quad tree. This is discussed for the case of 
use of ten bits of tile index. Each node in the quad tree can have only four children: none, some, or all of which may 
be set to nil. The node has pointers to these children. The four children (numbered 0, 1 , 2 and 3) correspond to spatial 
locations: top left, top right, bottom left and bottom right respectively. Each level in the tree therefore represents a 
spatial resolution double that of the previous higher level in the tree. In this case, there are ten levels in the tree - one 

so for each of the ten bits of tile index provided, separately, both in x and y dimensions. At each level, the appropriate bit 
from the y index indicates top or bottom (0/1 respectively) and the corresponding bit from the x index indicates left or 
right (0/1 respectively). In this way the 2 bit number formed by the y bit followed by the x bit gives the child node number 
at the current level in the Iree. "Leaf nodes" of the tree (those at the lowest level, level ten) are either empty (set to 
zero) or point to a tile directory which is represented as described above with reference to Figure 11 

ss Figures 12a and 12b show a 4 bit tree that represents the active and compressed tiles shown in figure 11 for the 
simple 2 dimensional tile index described above. Figure 1 2a is a spatial representation of the tile index which shows 
the different degree of spatial resolution that exist for different regions covered by the index. Only those regions which 
represent (or neighbour) scanned regions and have either active or compressed tiles are represented at the highest 
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resolution in the index. Figure 12(b) shows the index tree corresponding to the spatial index. Filled leaf nodes that 
point to tile directory entries are indicated by a vertical line beneath the corresponding leaf node in the tree. The tile 
entry 224 with 4 bit binary y index 0101 and 4 bit binary x index 1000 is shown as solid in Figure 1 2a: the corresponding 
array cell is shown as cell 204 in Figure 11 . The associated leaf node 225 is indicated in Figure 1 2b. As can be seen, 

s the binary numbers formed from each successive pair of bits from the two indices ( with the y bit in each case given 
first) give the node number at each level in the tree that lead to this leaf node. This is in the order 01 (nod e number 
1), then 10 (node number 2), then 00(node number 0) and finally 10(node number 1) again. 

Figure 13 shows a suitable representation of the tree as an array of quads 214: practical embodiments will have 
16 bits (i.e. 10 index bits, for both x and y dimensions), though for convenience the 4 bit case is illustrated in Figure 

io 1 3. The first quad in the array is always the root node of the tree. The locations within the array of all the other nodes 
in the current tree are dependent upon the order in which the tree was constructed. The construction of the tree will 
be described below. Before this, where we described indexing of the tree, it was assumed that the tree had been 
properly constructed and that the tile we wished to address was already in the tree. 

Figure 12 shows the bits used to index the tree arranged vertically with the y co-ordinate index on the left. As 

is discussed above, each pair of corresponding y and x co-ordinate bits forms a two bit address which corresponds to a 
child of a node in the tree. The first pair always represents a child of the root node (the first node in the quad array). 
The two bit address formed from the first pair of bits gives the number (0, 1 , 2, 3) of the 16 bit element within the quad 
that represents the child of the root node. The value stored at this 16 bit location gives the quad array index of the 
quad that represents the children of that node: there are again clearly four such children. The second pair of bits 

20 identifies the child at that level. This process continues in the same manner to the tenth bit pair (in the more limited 
case shown in Figure 13, the fourth bit pair) where the quad value now gives a tile directory entry number instead - 
this last node in the tree is a leaf node. 

The quad tree structure provides a very efficient representation of the tile index. Each leaf node of the tree repre- 
sents 4 tile directory entries. A tree of 700 quad nodes requires just 700*(4 children)*(2 bytes) to represent it; this is 

25 less than 6KB lor the whole tree. Most nodes in the tree are leaf nodes, as each set of four nodes has only one parent: 
there are % the number of parent nodes as leaf nodes; 1 /1 6 the number of grandparent nodes (in the sense ol parents 
of parents to leaf nodes); and 1/64 great grandparent nodes; and so on. These "ancestor nodes' can be summedas 
a series: the series converges to a number of total ancestor nodes which is less than 3/8 of the number of leaf nodes. 
In this case of 700 total quad nodes, we thus have more than 500 leaf nodes representing 2000 tiles: in the case of a 

30 scanner operating at 300 dpi, this is more than enough to scan the area of a whole letter sized page. A particular 
advantage provided by this data structure is that this area can be scanned in any direction. Furthermore, where larger 
arrays are used to represent the quad tree, the area that can be scanned increases linearly 

The construction of the quad tree is now described. This is not complex: all that is needed is an additional pointer, 
nextAvallableQuad, that points just beyond the highest quad array entry currently used. Initially only the root node is 

35 allocaled, all the pointers to its children are nil (represented by quad array index value 0) and nextAvallableQuad 
points at quad array entry 1 . When building the tree we begin to work our way down the tree as if we were indexing it. 
If at any level the two bit value composed fr om corresponding pairs of y index and x index bits indicates a nil pointer, 
then the current value of nextAvallableQuad is substituted and nextAvallableQuad is incremented to the next quad 
array location. The four entries for the children of the newly created node must of course be set to nil. In this way the 

w tree only grows during its construction, and no method to support deletion of nodes in the tree needs to be included. 
A method for tiling in the context of hand scanning is shown in the flow diagream of Figure 14. This includes the 
process of rectification, which occurs together with tiling: however, details of rectification processes are not discussed 
here, but rather at a later point in this description. End point position information is read by the tile manager 231 from 
the stream of increment data 230 ahead of the point at which it is used for rectification. The appropriate delay between 

45 the tile manager 231 and the rectification unit 237 is introduced by an increment hold buffer 232. This is used to hold, 
a predetermined number of increments that have been examined by the tile manager but have not been rectified. As 
will become clear from the discussion that follows, an appropriate number of increments to hold is (2" + 1 ), where n is 
an integer (typically 33 increments, in certain preferred embodiments). Once the increment hold buffer 232 is full, tile 
management is suspended unlil all the increments in the hold buffer 232 have been rectified. The last increment of 

so one such chunk of increments becomes the first increment of the next chunk (hence the combined tiling and rectification 
process advances in 32 increment steps for a 33 increment hold buffer 232). 

During each round of tile management the end point geometry associated with the increments in the increment 
hold buffer 232 (including the repeated increment from the previous round) is analysed to identify the set of currently 
active tiles. This can be done by considering the set of tiles through which a line drawn between the endpoints of each 

55 increment in the chunk passes. The superset of all such tiles (i.e. the combined group of all tiles intersected by an 
increment in the chunk) will be those that are active during the bout of rectification that will follow when the increment 
hold buffer 232 is complete. The set of tiles through which the line between the end point positions passes can be 
determined using a standard frame buffer line drawing technique: an example ol such a technique is discussed in 
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"Computer Graphics: Principles and Practice", James D. Foley, Andriesvan Dam, Steven K. Feiner and John F. Hughes, 
Second Edition, Addision-Wesley Publishing Company, Reading, Massachusetts, 1990, pages 72 to 81. Figure 15 
shows part of the space of possible tiles and the subset of active tiles that intersect the line drawn between end point 
locations of one increment. 

Before such a line can be drawn through the space of possible tiles, the end point position data must be scaled 
from rectilinear image pixel units to tile index units. In the case of 64x64 tiles as indicated to be appropriate for hand 
scanning, this involves division by 64. It is strongly desirable in practice to maintain all of the fractional precision that 
exists in the end point position data, in order to ensure that exactly the same set of tiles are made active as will be 
encountered during the forthcoming rectification process. 

It is assumed here that the space between two increments is not large enough to alio w a whole tile to lie between 
them. For the handscanning embodiment considered with use of 64x64 tiles this assumption will clearly hold, as the 
magnitude of the separation between increments is typically at most 2 output pixel units. It may in practice even be 
possible to sub-sample the set of increments (e.g. to consider for this purpose only every fourth one) used to identify 
active tiles with no loss of generality to achieve reduced computational overheads. 

As each active tile is made explicit in this way, it is added to a tile index table N of active tile indices. As a tile is 
encountered during the line drawing process, then the indices of that tile are compared against those present in the 
table N of active tile indices. If a compatible entry is already present, the new tile indices are ignored; otherwise they 
are added to the end of the table, and the table size is incremented by one. Each tile index pair requires 1 0 bits for the 
y index and 10 bits for the x index. In preferred embodiments, for reasons discussed further below, a pair of 16 bit 
unsigned integers is used to represent the tile index pair, with the y index again first and followed by the x index. 

The table N of currently active tile indices is used as the basis of a number of important tile management activities. 
These are: 

• identification of previously active tiles that are no longer active and can now be compressed; 

• updating the tile index 21 4 to ensure that new active tiles are present and have memory allocated for them; and 

• building a tile cache 236 to be used by the rectification unit 237 to give rapid access to the active tile data memory 
without need for reference to the tile index itself. 

In a particularly preferred embodiment of the method a pool of tiles is allocated in the Active Tile Data Memory 
Buffer 212. This consists of the two-dimensional arrays of byte data that represent the part of the image represented 
by each tile. The active tile data memory is partitioned into contiguous areas, each large enough to represent an image 
tile. The minimum number of active tiles required is equal to the largest number of tiles that can be indicated as active 
at any one time. This in turn is determined from the largest number of tiles that can lie beneath the chunk of increments 
in the increment hold buffer 232. Assuming that the spacing between increments is less than 2 output pixels, as should 
be the case for the handscanner embodiment considered, then the overall spacing of the chunk could be at most 32*2: 
64 output pixel units. When scanning at 45 degrees with respect to the rectilinear image the maximum horizontal extent 
will be 64*V2; 90 pixels. In a worst case, these span 3 tiles on each row The maximum number of tiles along the 
scanner is given by the length of the scanner in output pixel units L, divided by the tile size and truncated + 2 to allow 
for partial tiles (1 if L is an exact multiple of the tile size). This gives a conservative (as the scanner does not subtend 
its whole length at 45 degrees) maximum number of active tiles of 

maxActiveTiles = 3*(L/64 + 2) 

For a 2.5 inch scanner this gives 

maxActiveTiles2.5 = 3*(750/64 + 2) = 39 tiles 

for a 5 inch scanner it gives 



maxActiveTiles5 = 3*(1500/64 + 2) = 75 tiles 

and for an 11 inch scanner it gives 



maxActiveTiles1 1 = 90*(3300/64 + 2) = 159 tiles 
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assuming for each an output resolution of 300 dpi. Note that truncation in the division is assumed. 

If one has, as in the embodiment discussed here, a byte provided for grayscale levels, a 64x64 tile requires 4KB 
of memory, 2.5, 5 and 1 1 inch scanners at 300 dpi require a minimum of 156 KB, 300 KB and 636 KB of active tile data 
memory respectively (different levels of resolution, and consequently different memory requirements, are of course 

s entirely possible within the scope of the invention). These figures assume that all the previously active tiles that are 
no longer active are copied, compressed and stored in the Compressed Tile Data Buffer 213 prior to the creation of 
new tiles. As a result, the memory resources from the Active Tile Data Buffer 212 that represented tiles which are no 
longer intersected by the current chunk of increments are returned to the pool of active tiles. This allows the reallocation 
of active tile data memory to new active tiles via the tile index and associated tile directory entries. This approach 

10 requires that the compression and rectification processes work in sequence, with compression taking place first and 
followed by rectification. In practice, a parallel implementation is preferred. In a preferred embodiment the number of 
active tiles (and hence the size of the active tile data buffer 21 2) is increased with respect to the minimum requirement 
to allow a dual processor implementation. One processor and supporting hardware units performs rectification while 
another performs compression. The increase in number of active tiles required is determined by the minimum degree 

is of overlap that must be present between the successive chunks. 

We now describe how the Tile Manager 231 updates the tile index 214 to ensure that new active tiles have memory 
allocated to them. At the start of scanning the Tile Index 214 and Tile Directory 211 entries will indicate an empty quad 
tree. At this point, the array of quads that represents the Tile Index 214 has a single root entry at array location 0, all 
four children node pointers for the root entry are nil, and nextAvailableQuad points to quad array element 1 . An index 

to nextAvailableDirectoryEntry points to the first entry in the Tile Directory 211 . An Active Tile Data Offset FIFO 233 is 
initialised to include an entry for each block ol tile data in the Active Tile Data Buffer 212. When scanning begins, each 
entry in the new table of active tile indices N is examined, and the tile index associated with its location is accessed 
through the quad tree Tile Index 21 4. If the quad tree Tile Index requires a new entry to represent this tile, such a new 
entry is created as discussed above with reference to Figure 1 3. If an entry in the index table refers to a leaf node of 

25 the tree that is currently nil, this indicates that the scan path has entered a new tile. In this event, the leal node in 
question is made to reference the nexl entry of the Tile Directory 211 at the localion pointed to by nextAvailableDi- 
rectoryEntry: nextAvailableDirectoryEntry is then incremented. This new directory entry is marked to indicate that 
it has not been compressed, and it is made to point through its offset field to the Active Tile Data offset value extracted 
from the Active Tile Data Offset FIFO 233. The data associated with the tile image in the Active Tile Data Buffer 212 

30 is then initialised to a background value. This background value is chosen to be the brightest 8 bit white intensity 
available for the grayscale used. 

At this stage it can be determined which tiles were active in the last bout of scanning (in the case used in this 
description, last tranche of 33 increments) but are no longer active in 1he current bout. The indices ol such tiles are'in 
a preferred version of the method added to the Completed Tile Index FIFO 234. This approach requires extra storage 

35 within the Tile Manager 231 to retain the table of active tile indices for the previous chunk of 33 increments. At the start 
of the scan, this table will be empty. Each entry in the "stale" table of tile indices is compared against the newly acquired 
active table of indices for the current 33 increments of the Increment Hold Buffer 232. Any entries that are present'in 
the old table which are not also present in the new table are added to the Completed Tile Index FIFO 234 to indicate 
that they should be compressed by the Tile Compressor 236. 

io The Tile Compressor 238 extracts tile indices from the Completed Tile Index FIFO 234 in turn and compresses 
these tiles. Compression can include any appropriate form of processing to reduce the amount of memory consumed 
by a tile, including the toss of any component in the data whose retention is not considered necessary to achieve a 
representation of the desired quality or form. Compression could thus include changes of representation and/or spatial 
or greyscale resolution. An appropriate approach which is well known in the art is binarisation as a result of the appli- 
es cation of a thresholding scheme (such a scheme is discussed in "Digital Image Processing", William K. Pratt, John 
Wiley & Sons, Inc., New York, 1 991 , pages 597 to 600). Any one of a number of standard compression schemes for 
greyscale or binary data could be employed, such as group 4 fax compression for binary data or LZW for greyscale. 

During the compression stage, the original tile from the Active Tile Data Buffer 212 remains unaffected. A new 
copy of the tile in compressed (or otherwise changed) form Is made within a Compressed Tile Data Buffer 21 3 location 

50 indicated by a nextAvailableCompressedData pointer. This buffer is filled sequentially in the order that tiles are 
marked for compression. The Tile Directory 21 1 entry associated with the now-compressed tile (obtained from the Tile 
Index 235) is then updated to point at the appropriate location within the Compressed Tile Data Buffer 213 and both 
its length field and compression state are amended. The offset within the Active Tile Data Buffer 212 associated with 
the active version of the now-compressed tile is then returned to the back of the Active Tile Data Offset FIFO 233 to 

5£ be used when required to represent a further active tile. Finally, the nextAvailableCompressedData pointer is incre- 
mented in accordance with the length in bytes of the compressed tile. In general, each tile will compress to a different 
number of bytes, depending upon the content of the tile and the compression scheme employed. 

As described above, a Tile Cache 236 is provided for computational efficiency in the rectification process. Figure 
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1 6 shows the Tile Cache built by the Tile Manager 231 . To address a specific pixel at location (x, y), where y and x are 
each 16 bit row and column positions, each co-ordinate must first be broken into its tile index component, and its "within 
tile' position component. These are provided by the top 10 bits and the bottom 6 bits of each co-ordinate respectively. 
The way the y and x 10 bit tile indices are used to access the tile array data within the Active Tile Data Memory Buffer 
s 212 using the Tile Cache is shown in Figure 16. The y index offset, corresponding to the minimum y index amongst 
the currently active tiles, is subtracted from the y index o( the co-ordinate under consideration to give a local y index 
that is used to access both an x index offset table 241 and a tile data offset table 242. The x index offset table 241 is 
a linear table of x tile offsets which contains the minimum x index for each row (and hence y tile index value) of the tile 
index after the y offset row, up to a maximum row determined by the length of the scanner (this "max row" will be y 
10 index offset + (L/64) + 2, where once again the added factor of 2 results from the division and the consequent effect 
of partially filled tiles and truncation). The tile data offset table 242 is an array of triples, each element of which contains 
a 1 6 bit index to the Active Tile Memory Buffer 2 1 2. The entry in the x index offset table 241 is given in the case shown 
by adding the local y index 1o a table base address (assuming that arithmetic is performed in x index offset table cell 
units: in this implementation, this will be 16 bit). Similarly, the relevant triple in the tile data offset table 242 is given by 
'5 adding three times the local y index to the tile data offset table's base address. The index into the tile data offset table 
242 is further altered so as to address the relevant tile index by adding the local x offset for the relevant row: this is 
given by Ihe difference between the original x tile index and the entry for that row from the x index offset table 241 . 

The content of the relevant cell in the tile data offset table 242 provides a pointer to the base address of the array 
of image data that represents the tile within the Active Tile Data Memory Buffer 212. The low 6 bits from each co- 
ze ordinate are then used to form an address within the tile in the usual way. Not all cells in the tile data offset table 242 
will be valid: those that are not have anil pointer. The table is constructed to be large enough for the worst case scenario 
(in terms of tiles required by a chunk of, here, 33 increments) and frequently fewer tiles than can be indexed by the 
tile data offset table 242 will be currently active. 

Rectification processes suitable for use with the combined tiling and rectifying scheme discussed above with ref- 
2S erence to Figure 14 are discussed below. It would be possible to adopt a Bresenham algorithm approach to achieve 
rectification, as is discussed in International Patent Application Publication No. WO 96/27257 and US Patent No. 
5578813. However, more sophisticated approaches such as those discussed in European Patent Application No. 
97302519.0, the content of which is incorporated within this application to the extent permitted by national law, are 
preferred. 

30 Figure 1 7 shows a coordinate frame provided by the array of tiles to wh ich the stream of increments 88 are mapped 
to form the final rectilinear image. There is no rectilinear image buffer of predetermined dimension, as in the methods 
of International Palent Application Publication No. WO 96/27257 and US Patent No. 557881 3, but rather a tiled recti- 
linear image space of undetermined overall dimension but predetermined scale. The X and Y co-ordinates of the 
position tags associated with each increment correspond to horizontal and vertical dimensions (column and row indices) 

35 in the rectilinear image space. Also shown are the loci of positions tags (shown enlarged in the scrap view for a short 
section, on the right) for each end point of the linear image sensor associated with a stream of increments. Within 
these a pair of increments (Increment #1 and Increment #2) are highlighted by joining their position tags with straight 
lines. These are chosen as almos t intersecting in the overlap region where the scan doubles back on itsell. 

It is possible that the resolution of the rectilinear image buffer, which is determined by the output resolution required 

40 from the scanner (typically either 200 or 300 dpi) differs to that in which the endpoint position data is measured. That 
is in turn determined by the resolution of the navigation sensors 24 (which is dependent upon the spatial extent of the ■ 
features, e.g. paper fibres, imaged during navigation). To accommodate such differences it is necessary for the endpoint 
position data to be scaled to the output pixel resolution. 

The next operation is to map the position-tagged image increments within a swath into the tiled rectilinear image 

45 space. It is particularly important that sufficient image quality is maintained in this process. One approach is simply to 
map the elements within each increment to the pixels of the rectilinear image array through which the straight line 
drawn between the endpoint locations passes. The mapping could either be to the closest linear sensor element for 
each image pixel, or involve interpolation between sensor elements. Appropriate linear or cubic interpolation methods 
for application in this context are described in Wolberg, "Digital Image Warping", pp 127-131 , IEEE Computer Society 

so Press, Los Alamitos, CA, 1 992. 

Performing rectification in this way will inevitably introduce sampling artefacts because each increment has a 
distinct impact. Moreover, it will be necessary in practice that the original swath of position tagged image data is strictly 
over-sampled with respect to the rectilinear pixel grid (which may, for convenience, be at a higher spatial resolution 
than the pixel pitch of the linear sensor itself) otherwise pixel drop out will occur in the rectilinear image. A solution is 

55 illustrated in Figure 1 8. In this ligure, the end point positions of a pair of increments 11 and 12 define a region. These 
increments can be a successive pair of linear sensor readings - the use of consecutive readings is termed here the 
PAIRWISE approach. Alternatively the increments 11 and 12 can be the bounding pair of a larger group of sensor 
increments - the use of such increments is here termed the CHUNKWISE approach. It will be shown that there is 
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computational advantage to be gained in the CHUNKWISE approach from treating a group of increments in unison at 
the cost of some small reduction in image quality when compared with the PAIRWISE approach. 

Pixel locations in the tiled rectilinear image space are associated with integer locations. For each such pixel within 
the region defined by a pair of increments, an intensity value is then computed by integrating a number of linear sensor 
s pixels that map to a region around it. In particular preferred embodiments, various forms of bilinear interpolation are 
used to integrate a specific four linear sensor pixels taken from the pair of successive increments that span the pixels 
location: two sensor pixels from the increment that maps to a line that lies to one side of the rectilinear array pixel, and 
a further two sensor pixels from the adjacent increment that maps to a line on the opposing side of the rectilinear array 
pixel. 

10 In Figure 18 a pair of increments 11 and 12 are shown mapped into a portion of the rectilinear space that defines 
the scanned image. For simplicity end point positions P1 , P2 and PI ', P2' are made to correspond to the centre of the 
first and last sensor elements of increments 11 and 12 respectively. Intermediate sensor elements such as Al, A2 and 
B1 , B2 are also shown. Pixels locations in the rectilinear image array are chosen to correspond to integer locations 
with respect to the end point position geometry. The resolution of the output rectilinear image may not be the same as 

is that of the linear image sensor. However, irrespective of the actual resolution, N sensor pixels map along the line that 
joins the endpoints of each increment. For ease of explanation, in Figure 1 B N=8, that is each linear sensor has just 8 
elements in it and the length of each increment is N-1 =7 times the pixel pitch of the linear sensor as measured in output 
rectilinear pixel resolution units. 

Specific pixels that lie within the region defined by a bounding pair of increments can be identified from the interval 

20 given by the x components of the intersects of each rectilinear image row with the two increments. An example is 
provided by pixels C1 and C2 on rectilinear image row D in Figure 1B. From simple geometry, if increment 11 has 
endpoint co-ordinates (x1 , y1) and (x2, y2) then the x co-ordinate of its intersect with row D which has integer value 
yD is given by 

xD = (yD-y1 )*(x2-x1 )/(y2-y1 ) 
it follows that the intersection with row D + 1 is given by 

30 

x(D+1) = xD + (X2-X1 V(y2-y1) = xD + xlnc 

where xlnc is a constant increment for each successive row. Similarly if increment 12 has endpoint co-ordinates (xV, 
y1 ') and (x2', y2') then the x co-ordinate of its intersect with row D which has integer value yD is given by 

3S 

xD' = (yD-yl'nxa'-xI'HyE'-yV) " 
it follows that for the intersection of row D + 1 

40 

x(D+1)' = xD" + (x2'-x1 ')/(y2'-yV) = xD' + xlnc' 

This provides a computationally effective way to compute the range of intersection for each row in turn. The intersections 
45 with the current row are computed by simply incrementing the intersections with the previous row. 

First we must specify the range of rows that lie within the region. This is given by the range of integer y values that 
are common to both increments. These are the integer values in the range 

so [Ceiling(Max(y1 , y1 ')) , Floor(Min (y2, y2'))] = [Y1 , Y2] 

assuming that the orientation of the scanner with respect to Ihe paper remains roughly vertical. This approach can 
easily be extended to allow the scanner to be at an arbitrary orientation with respect to the page. For example, if the 
scanner is inverted the linear sensor data can be reversed and the polarity of the end point position data reversed, 
ss Furthermore, if the angle of the scanner is more than 45% from vertical correct rectification can be achieved by flipping 
the polarity of both the x/y position co-ordinates and the rows and columns of the image. For ease of explanation, the 
discussion of the embodiments will only deal with the roughly vertical case, but extension to the more general cases 
indicated above is a straightforward matter for the person skilled in the art. 
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pixels in the region can be identified in turn using the following pseudo code 



TRANSFORM LOOP 1 
{ 

// INITIALISATION 
y = Yl 

xlnc = (x2-xl)/(y2-yl) 
xlnc' = (x2'-xl')/(y2'-yl>) 
xD = xl + (y-yl)*xlnc 
xD' = Xl' + fy-yl^+xlnc 1 

// MAIN LOOP 
while (y<=Y2) 
{ 

El, X2] = [Ceiling(xD) 5 Floor(xD')] 
x = X1 

//INNER LOOP 
while (x<=X2) 

VISIT PIXEL (x, y) 

xD + = xlnc 
xD' + = xlnc' 

y++; 

} 



whe re operators floor and ceiling have their usual mathematical interpretations. The actual details of the initialisation 
of^ 

thil/n thf ! ' h6 assum P ,,on 13 made ,hat xD ^ always less than xD'. It is straightforward to extend 

this to the generic case by testing the order of xD and xD' on each row. 

^J h lZ« Sl6P i S t0 VT he WhiCh ' inear S8nSOr piX6ls make a con,ribu,i °" 10 rectilinear grid pixel and in 
what proportions. An embodiment of the PAIRWISE approach is now described 

In this embodiment of the PAIRWISE approach to the method according to the invention, recreation is based 

e" zzsnz: ? in rr er 9raphics texture mappins (see ,or ^ ^z^t 

fncremen^ 1 Z ?t ™ " 1 T 'T M d, " mitod ^ ^ tW ° lineS ,hat j° in ,he end P° int Prions of the 
increments 11 and 12 ,s mapped backwards into a sub pixel location within a coordinate space defined between the 

«nZf T^ 65 - ThiS Si,Ua,i ° n ' HIUStra,9d fi9Ure 1 9 the p - <* "ncr-mSTlfaSSSS 
spacing, the increments align wrth the y axis and the increments are placed a unit pixel apart Wrthin this space the 
top sensor element of the left hand increment has co-ordinates (0 0) P 

ordiTx 6 ^ TZTv TT ima9e " Sh ° Wn maPP6d be,Wee " the incremen,s 11 and 12 11 has ^b-pixel co- 
n ST ♦ a ,hB ' n1e9er C0mp0nent ° f ,he C0K3rdinate alon 9 increment and x and y areThe sub- 

i s That ^. r6preSenta,n/e PiXe ' V3lUe C9n then be d *<— -ing bilinear interpolation o, the four s roundlg 



BiLinear(M ,l2,Y,ypc) = H[Y]*(1-x)*(1-y) + I2(Y]*(1 -y)*x 
+ H[Y + ir(1-x)*y + l2[Y + 1]«x*y 
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The best approximation to the transform from the rectilinear image to the space between increments involves the 
solution of a quadratic for every pixel. This is nol a satisfactory practical solution because of the expense of computation 
involved. While it is possible to derive approximate linear homogeneous mappings from rectilinear image coordinates 
to the co-ordinate space defined by a pair of increments, a non-homogeneous solution is preferable. With such a non- 
homogeneous solution no assumptions about the nature of the local transformation are required. It is also possible to 
solve easily for corner cases, such as when the projection of an adjacent pair of increments involves a cross-over in 
the rectilinear image space. 

In such a non-homogeneous method, a different mapping is defined for each row of the rectilinear image in the 
interval [Y1 , Y2], This is again determined from the intersection of that row with the lines that connect the end-point 
positions. In this case, the distance of the intersection point along the line itself is recovered. This corresponds, after 
any scaling required for the differences in resolution, to a physical position along the linear sensor. In Figure 19 aD 
and aD' are the intersection along increments 11 and 12 respectively for rectilinear row D; where from simple geometry 

aD = (yD-yir(N-iy(y2-y1) 

and 

aD' = (yD-ylT(N-iy(y2'-yV) 

it follows for the intersection of row D 4 1 

a(D 4 1 ) = aD + (N-1 )/(y2-y1) = aD 4 aRowInc 

and 

a(D4l)' = aD' 4 (N-iy(y2'-y1 ') = aD + aRowInc' 

which again gives a simple recurrence relationship for each successive row for efficient serial implementation. 

Each pixel in the rectilinear image along row D in the interval [X1 , X2] maps to a different point along a line in the 
increment space that joints the two points of intersection at (0, aD) and (1 , aD') respectively. Assuming a linear mapping 
along this row the pixel at location (xE, yD) in the rectilinear image will map to location (aE, bE) in the increment space 
where: 

aE = aD 4 (xE-xD)*(aD'-aD)/(xD'-xD) 

and 

bE = (xE-xD)/ (xD'-xD) 
and again it follows lor the next pixel along row D that 

a(E4l ) = aE 4 (aD'-aD^xD'-xD) = aE 4 alnc 

and 

b(E4l ) = bE 4 1/(xD'-xD) = bE 4 bine 

Where the intensity at (aE, bE) is determined by the bilinear interpolation. 

The additional initialisation and modified inner loop of TRANSFORM LOOP are 
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TRANSFORM LOOP 2 
{ 

// INITIALISATION 
y = Yl 

xlnc = (x2-xl)/(y2-yl) 
xlttc' = (x2'-xl')/(y2'-yr) 
xD = x1 + (y-yl)*xlnc 
xD' = xl'-Kj-yl'^xInc' 

aRowInc = ( N ■ 1 ) / ( y 2 • y 1 ) 

aRowInc' = (N-l)/(y2'-yl') 

aD = (y-yl)*aRowInc 

aD' = (y-yl')*aRowInc' 



// MAIN LOOP 
while (y<=Y2) 
{ 

[XI, X2] = [Ceiling(xD), Floor(xD')] 
x = Xl 

alnc = (aD'-aD)/(xD'-xD) 
bine = l/( xD'-xD) 
a = (xl-xD)*aInc 
b = (xl-xD)*bInc 
//INNER LOOP 
while (x<=X2) 
{ 

A = Floor(a) 

pixel[y][x] = BiLinearfll, 12, A, a-A, b) 
a + = alnc 
b + = bine 

} 

xD += xlnc 
xD' + = xlnc' 
aD 4- = aRowInc 
aD' + = aRowIInc' 

y+ + ; 

} 

} 



More direct geometrical interpretation can also be used in other embodiments of the method according to the 
invention to perform PAIRWISE rectilinear pixel interpolation. Such methods do not require the explicit mapping of 
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rectilinear image pixels back into an increment co-ordinate space. 

One particularly simple embodiment having reduced computational requirements with respect to the back projec- 
tion approach is to perform interpolation along each increment and then interpolate these values directly in Ihe rectilinear 
image. Better image quality can be achieved, at the cost of increased computation, by mapping each rectilinear image 

s pixel to the closest point on each of the increments. These are the points on the increments whose interpolated (along 
the linear sensor) intensity best reflects the intensity of the rectilinear pixel under consideration. The pair of interpolated 
increment intensity values are then combined, again using a linear interpolation, to give the best representative intensity 
for the pixel. Such embodiments are described in European Patent Application No. 9730251 9.0. 

Next further embodiments of a rectification approach employing the CHUNKWISE method of rectilinear mapping 

10 will be presented. The CHUNKWISE approach has the twin advantages that the outer loop of the appropriate version 
of TRANSFORM LOOP is computed less frequently and that larger numbers of output pixels are accessed sequentially 
which leads directly (without the need for a cache) to improved memory bandwidth. 

Figure 20 shows a chunk consisting of M, in this case 4, increments. The outer pair of increments, labelled 1 1 and 
12, are again used by the outer loop of TRANSFORM LOOP to identify valid pixels that map within the chunk. Preferred 

'5 embodiments determine rectilinear array pixel values with an extended form of non-homogeneous back projection 
algorithm as described previously lor the PAIRWISE case. Figure 21 shows the transformation of a row labelled D in 
Figure 20 into an increment space that contains the 4 increments bounded by 1 1 and 12. It joins the points of intersection 
of row D with each of the outer increments. Following on from the PAIRWISE case, as shown in figure 19, these 
intersections have offset distances along the linear sensor of aD and aD' respectively for 11 and 12. The co-ordinates 

20 of the intersection points in the increment space are thus (0, aD) and ((M-1), aD'). 

Each pixel in the rectilinear image along row D in the interval [X1 , X2] maps to a different point along a line in the 
increment space that joins the two points of intersection at (0, aD) and ((M-1), aD') respectively. Assuming a linear 
mapping along this row the pixel at location (xE, yD) in the rectilinear image will map to location (aE, bE) in the increment 
space where: 

25 

aE = aD + (xE-xD)*(aD'-aD5/(xD'-xD) 



bE = (M-1)*(xE-xD)/(xD'-xD) 
and again it follows for the next pixel along row D that 

a(E+1) = aE + (aD'-aDV(xD'-xD) = aE + alnc 



b(E+1 ) = bE + (M-1 y(xD'-xD) = bE + bine 



Where the intensity at (aE, bE) is determined by the bilinear interpolation. 
45 Hence TRANSFORM LOOP becomes: 



55 
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TRANSFORM LOOP 3 
{ 

// INITIALISATION 
y = Yl 

xlnc = (x2-xl)/(y2-yl) 
xlnc' = (x2'-xr)/(y2'-yr) 
xD = xl+(y-yl)*xlnc 
xD' = xl' + (y-yl')*xlnc' 

aRowInc = ( N ■ 1 ) / ( y 2 ■ y 1 ) 

aRowInc' = (N-l)/(y2'-yl') 

aD = (y-yl)*aRowInc 

aD ' = (y-yl')*aRowInc' 



// MAIN LOOP 
whfle(y<=Y2) 
{ 

[XI, X2] = [Ceiling(xD), Floor(xD')] 
x = Xl 

alnc = (aD'-aD)/(xD'-xD) 
bine = (M-l)/( xD'-xD) 



a = (xl-xD)*aInc 
b = (xl-xD)*bInc 
// INNER LOOP 
while (x<=X2) 
{ 

A = Floor(a) 
B = Floor(b) 

pixel[y][x] = BiLineardB, IB+1, A, a-A, b-B) 
a + = alnc 
b + = bine 



xD + = xlnc 
xD' +=xlnc' 
aD + = aRowInc 
aD' + = aRowIInc' 

y++; 
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Where IB is the Bth increment in the chunk and IB + 1 is the (B + 1)th increment in the chunk. 

The CHUNKWISE method relies on the fact that scan path is almost uniform and only smoothly changing. This 
must be true both spatially (the scan should be a uniform shape) and in terms ol the sampling frequency along the 
scan path (the scan should be sampled uniformly). The CHUNKWISE method is particularly well adapted for use with 

s the tiling method illustrated in Figure 1 4. If the increment hold buffer 232 is adapted to hold (S + 1 ) increments, then 
uniform chunks with (T + 1 ) increments can be employed, where T is an integerfactor of S. For example, if 33 increments 
are held in increment hold buffer 232, then uniform rectification chunks can be formed of any of 2 (the PAIRWISE limit), 
3, 5, 9, 17 or 33 increments. The rectification step will then occur S/T times between each successive replenishment 
of the increment hold buffer 232. 

10 Also possible is the extension of the CHUNKWISE method by variation of the chunk-size dynamically depending 
upon the degree of uniformity of the scan path. Where uniformity is good, a large chunk can be accommodated without 
introduction of image distortion, but where the scan path or sampling regime were to vary rapidly, only a small chunk- 
size could be accommodated without introducing artefacts. In the limiting case, such embodiments revert to a PAIR- 
WISE method. 

'5 As has been discussed above, the rectilinear image space available to memory can easily accommodate a single 
letter sized or A4 sized page at the resolution of the linear sensor (typically either 200 or 300 dpi). As there is no 
dimensional restriction to the rectilinear image space, the final image will be well formed independent of the initial 
starting position or the scan path. To obtain proper alignment, either the scan must begin at the assumed orientation 
(e.g. always parallel to the side of the page) or orientation must be recovered from the content of the scan and be used 

20 to reorient the final image. Methods for automatically determining the dominant orientation of text on a pagB, which 
can be used as a basis for the latter, are known in the literature eg. "The Skew Angle of Printed Documents" by H.S 
Baird. Proc 4th SPSE Conference Symposium on Hybrid Image Systems, Rochester, New York 1 987. 

The next step is to stitch successive image swaths within their region of overlap. The aim is to combine multiple 
swaths so that the result is a final rectified reconstructed image stored in tiled form - advantageously, this image will 

25 also be compressed. This should be dona in such a way as to identify and correct most of the accumula ted navigation 
error and to m-ask any residual error. 

In the embodiments to be described, the stream of navigation data provides the registration information needed 
for stitching. Since the navigation signal tends to accumulate error, it is continually amended by feeding back a cor- 
rection signal derived from analysis of feature offsets. 

30 However, first we will describe a method of stitching image swaths which assumes that there are no navigation 
errors. 

Some area of overlap is necessary in order to stitch two image swaths. Swaths are delimited by a reversal in the 
path of the scanning device back over a region of the original of which a part has just been scanned. A swath comprises 
the image data captured during a scanning swipe across an original being scanned. In the folbwing description, the 

3S term 'swath' will also sometimes be used to refer to the part of the reconstructed image formed from mapping such data. 

The scan represented by the loci of endpoints of increments in Figure 17 generates two overlapping swaths. In 
Figure 22 the portion mapped into the rectilinear image space that comprises Swath #1 is partially re-mapped on the 
return by the portion of the scan that corresponds to Swath #2, with the area of overlap 102 shown filled. At a time T 
a partial swath has thus far been scanned. Satisfactory stitching can be achieved, in this case, by simply continuing 

*o to map Swath #2 1o the rectilinear image space on top of Swath #1. Returning to Figure 17, Increment #2 can be 
mapped along its whole length in just the same way as for Increment #1 and similarly for each pixel in the region of 
overlap 1 02 in Figure 22. The final rendition in the rectilinear image space that resulted from Swath #1 is simply replaced 
with that resulting from Swath #2. 

In this event, the method set out with regard to Figure 14 will need to be amended to allow the content of a no 

45 longer active, and hence compressed, tile to be replaced - this is likely to require additional memory for the storage of 
compressed tiles, as if memory is allocated to compressed tiles by the tile compressor 238 according to the memory 
required by each tile, difficulties will result if the second version of the tile requires more memory after compression 
than the first version. Difficulties may also occur in the case of tiles only partially filled in the first swath, but rendered 
inactive (tiles at the edge of the linear sensor), ff these are compressed on the first pass, a second "full" version will 

so clearly take up more memory. A possible approach is never to compress a tile which is not filled - partially filled tiles 
can be checked for by evaluation of each corner of the tile to determine whether information has been written to it at 
the compression stage. Partially filled tiles could then be deleted when they become inactive (relying on swath overlap 
to obtain a complete image), or stored in a separate memory area. 

A scheme in which pixels in the area of overlap 102, mapped as filled tiles in Swath #1, are not replaced by those 

55 mapped by Swath #2, will pose fewer difficulties in the tiling approach, provided that all pixels mapped by Swath #2 
that are not part of the overlap (i.e. are nol also mapped by Swath #1 ) are properly mapped. That is to say, the region 
mapped by Swath #1 is used to clip the region mapped by Swath #2. In principle, provided that Swath #1 and Swath 
#2 correctly map all pixels that lie outside the overlap region, then the pixels in the overlap region can be obtained 
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squally from either Swath #1 or Swath #2, or a combination. 

In practice, due to the accumulation of errors between the position tags of Swath #1 and Swath #2, this simple 
approach which assumes no navigation errors does not give very good results. 

An embodiment of an approach to stitching in accordance with the present invention will now be described with 
reference to Figs 23 and 24. Figure 23 indicates the processing stages and data buffers involved while Figure 24 relates 
to the process utilised with respect to Swath #1 and Swath #2. The image increment data is mapped to the tiled 
rectilinear image space as described previously. Figure 24 shows use of the image captured in Swath #1 in mapping 
on to the rectilinear image space of the return pass, Swath #2. Navigation correction is calculated by correlating features 
within the area of overlap between Swath #1 and Swath #2. 

Fig. 24 highlights thisoverlap area which is divided into two parts 104 and 1 05. As shown in Fig. 24, during collection 
of Swath #1, quadrilateral image segments 130 (henceforth called "registration tiles") are periodically labeled along 
the lower edge of the swath in the overlap area 105. On a later pass (Swath #2) the overlap area 104 of Swath #2 
above the area 105 containing the registration tiles 1 30 of Swath #1 is clipped, i.e. discarded, as Swath #2 is acquired, 
by not permitting those filled tiles provided by Swath #1 to be overwritten. The registration tiles 1 30 from Swath #1 are 
located in the top of what remains of Swath #2, after the overlap area 104 has been clipped. If the navigation data is 
perfect, there will be no offset between the location of a registration tile 1 30 and the location of that tile's rescanned 
image in Swath #2. More realistically, some navigation error will have accumulated since the last registration was 
performed. The offset between the appearance of this tile in the two swaths produces a correction factor which is then 
used to update future navigation position-tags associated with the image data, in order to minimize the total accumu- 
lated error. In this way the total accumulated error in the navigation data is prevented from growing so large that it 
introduces an obvious distortion in the region where the swaths overlap. 

Navigation error can be minimised by beginning to stitch Swath #2 to Swath #1 from the point where these swathes 
are closest in the stream of increments. The navigation errors should be least between these closely spaced increments 
in a scan sample. Continued use of matching around a turning point" between successive swathes (where obtained, 
as is advantageous, according to the capture approach shown in Figure 1 ), together with progressive use of detected 
error between swathes being stitched together to predict and process error further along the swathes, allow sell-guiding 
control of the stitching in a manner essentially analagous to the closure ol a zip or zipper. 

The processing stages for stitching Swath #1 to Swath #2 will now be described with reference to Figs 23 and 24. 
Figure 23 shows the image increment buffer B and the tiled rectilinear image space 100. Fig 23 also shows a feature 
location buffer 1 31 and an image buffer 1 32. An arrangement o* this type is shown in International Patent Application 
Publication No. WO 96/27257. However, in the tiled scheme described here an additional feature data buffer 135 is 
required. When a registration tile 1 30 is constructed, it is necessary to obtain a copy of the image tile containing that 
registration tile. This is because when such image tiles become inactive, they will be compressed by compresser 23B 
and the data they contain will no longer be readily available for stitching. It is thus advantageous to copy these tiles in 
uncompressed form to the additional feature tile buffer, with a pointer to each relevant tile in the additional feature data 
buffer 1 35 from the feature location buffer 1 31 . The processing stages are as follows: 

1. As mentioned above, during collection of Swath #1, registration tiles (130) are periodically labeled along the 
lower edge of the swath in the overlap area 105. An entire registration tile could be used for the correlation decribed 
above, but in the preferred embodiment a small area of high-frequency contrast (henceforth called a "registration 
feature") consisting of a rectangular tile (e.g., 15x15 pixels) of grey scale image is located within a registration tile 
forming part of the image being reconstructed in the tiled rectilinear image space 100 as a result of the capture of 
Swath #1 . Image tiles containing these registration tiles 1 30 are copied to the additional feature data buffer 1 35 before 
such image tiles are compressed, in embodiments where compression is provided as indicated in the method of Figure 

2. The position tags of the registration features (which define the position of each registration feature within the 
tiled rectilinear image space 100) are saved in the feature location buffer 131 prior to the start of the mapping of Swath 
#2 into the rectilinear image space 100. 

3. As Swath #2 is mapped, registration feature locations are identified in advance of that region of the rectilinear 
image space being written to by Swath #2. This is achieved by defining a rectangular capture window 107, indicated 
in Fig. 24 with a length equal to an image increment and a width of a few pixels, which precedes the image increment 
of Swath #2 which is currently undergoing mapping into the rectilinear image space 100. When a registration feature 
location stored in the feature location buffer 1 31 falls within the capture window 1 07, the registration feature location 
is selected (only one registration feature location may be selected at any one time). 

4. The selected registration feature location is used to define the predicted location of the registration feature when 
Swath #2 is mapped into the rectilinear image space. Navigation errors found with respect to previous registration 
features are taken account of by using the current error estimate (stored in an error buffer 1 34) to predict the location 
of the current registration feature. Hence, the search area used to locate the current registration feature then need 
only be large enough to account for the possible error increment. 
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5. The special image buffer 1 32 is used to store temporarily, directly from increment buffer B, a rectinlinear image 
mapping of the image data of Swath #2 located around the selected registration feature location as well as storing the 
registration feature location. In other words, at time T, the image increment 13B shown in Fig 24 is mapped into the 
buffer 1 32. The size of the special image buffer 1 32 must be adequate to store the registration feature plus the required 

5 search area, the size of which is calculated according to known error estimation algorithms. 

6. The image fragment from Swath #2 stored in the special image buffer 132 is then compared with the corre- 
sponding image fragment stored in the feature data buffer 135 pointed to by the feature location buffer 131 for the 
relevant registration tile 1 30. In this way, an offset is obtained which, when accumulated with the previous error estimate, 
gives an updated estimate of navigation error. This updated estimate is stored in the error buffer 134 along with a 

10 position tag. 

The position tag of the error estimate is simply the location, with respect to the rectilinear image space, of the 
centre of the current registration feature. It is used in turn to determine the increment in the increment buffer B to which 
the error estimate relates (i.e. the first increment corresponding to Swath #2 which maps to that image location). This 
indicates the increment by which the measured error estimate should be fully accommodated. 
is While other correlation approaches could be employed, an acceptable approach for calculation of the offset be- 
tween two image fragments is a "sum of squared difference" correlation. A small search area is defined around the 
original location of the feature and correlation coefficients are determined by equation: 

zo c k) = ' '»k,jti) 

where Ty denotes the grey scale values of the feature from Swath #1 and l i+kj+ | denotes the grayscale values of 
the newly acquired feature from Swath #2. Indices i and j specify pixel locations within the features, while k and I specify 
the magnitude of the proposed translational offset (constrained to remain within the search space). The smallest ele- 

25 ment in the resulting correlation array denotes the offset between the two features. The registration feature within the 
registration tile is selected to maximize image variance, since this improves the accuracy of the correlation method. 
Details of the application of this approach are discussed further in International Patent Publication No. WO 96/27257. ... 

In a second processing phase, the image data from Swath #2 is written into the rectilinear image space taking ? 
account of the error estimates at the recorded positions so as to remove those errors from the final reconstructed 

30 image. This is done by modifying the position data associated with the endpoints of the individual image increments 

of sensor data. This processing is advantageously completed, in the context of the Figure 1 4 method, before deter- ' 
mining which tiles are required as active tiles in a "bout" of tiling, rectification and compression. 

This preferred approach gives improved stitching because errors in the positions of the navigation sensor travelling 
through the overlap region between swaths are identified and corrected for before the final image is reconstructed. In * 

35 addition, error between the calculated positions of the top and bottom of the image increments which has accumulated ? 
during collection of the previous swath can be absorbed at once at the beginning of the next swath without introducing 
unwanted artefacts. This is indicated in Fig. 24 by the discontinuity between the right hand edges of the turning point 
between Swath #1 and Swath #2 derived from the first and second processing phases described above. 

The stitching step may be implemented in different ways with regard to the nature of image capture (more specif- 

40 ically, but not exclusively, scanning) required. One alternative is to require scanning to be carried out from the top of 
a page to the bottom, in which case stitching need only ever be done between the bottom of one swath and the top of 
the next. Another approach is to allow scanning to start at any part of the page but to require the initial direction of 
scanning to be maintained. In that case, there must be the capability to identify registration features on both edges of 
swaths, but once the direction of scanning is established, error estimates need only be kept for one side of the current 

*s swath. In another approach, scanning may be permitted in any direction and changes in scanning direction can also 
be accommodated eg. a spiral scan. In this third approach, not only must there be the capability to identify registration 
features on both edges of swaths, but error estimates must be kept for both the top and bottom of each swath in case 
the direction of scanning should change. The latter approach gives the greatest flexibility for the user but has a higher 
computing overhead. 

so In a preferred embodiment, the processing electronics for image reconstruction, stitching and image management 
is contained within the housing that defines the scanning device 10 of Fig. 1 . Thus, the scanned image may be imme- 
diately presented at the image display 16. However, the scanning device may contain memory to store the position- 
tagged image data, but without processing and file management electronics and firmware. 

As noted in reference to Fig. 3, the navigation and imaging sensors 22, 24 and 26 are preferably mounted on a 

ss pivoting member 20. In one embodiment, the pivoting member is connected to the remainder of the housing by at least 
one elastomer for which one end of the elastomer is connected to the stationary portion of the housing and the other 
end is connected to the pivoting member. The elastomer acts as a hinge. Thus, the pivoting portion is allowed to float" 
without the use ol frictional elements. Power, control and data signals may be conducted to the sensors via flex cables 
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that are shielded in order to minimize electromagnetic interference. Other methods of pivotalV attaching the pivoting 
member can be used. If the pivoting member is deleted and the sensors are in a fixed position on the housing care 
must be taken no to tilt the scanning device 10 excessrvely during image capture. In this embodiment, the design of 
illumination and optical elements must be given increased attention. 

While the invention has been described and illustrated as one in which a planar original is scanned, this is not 
critical. In fact, persons skilled in the art will readily understand how many of the techniques may be used for scanning 
three-dimensional images . However, the preferred embodiment is one in which the image of interest is formed on a 
medium, such as a piece of paper, a transparency, or a photograph, and the scanning device b in contact with the 
medium. 

As indicated above, however, in another aspect the invention is equally applicable to image capture of a different 
kind. The invention can equalfy be applied to capture of image data with an area sensor, such as a conventional digital 
camera. For capture of a large image at increased resolution with a digital camera, an effective approach is to scan 
across the image taking a series of image frames each capturing a part of the large image, oversampling so that there 
will be areas of overlap between different image frames captured by the digital camera, and then establishing relative 
positions between the different elements by using these areas of overlap. From these relative positions the image can 
then be reconstructed according to the method of the invention - no separate position sensor is thus required It would 
of course be possible to provide a position sensor with an area image sensor, but in this case the first aspect of the 
invention described wrih respect to the preferred embodiment is entirely applicable - however, in some capture regimes 
use of such a position sensor may be expensive, inconvenient, or even impossible 

The step of determining relative positions will now be discussed. While the skilled man will appreciate that it is 
possible to address more complex situations, the most straightforward applications of this approach are to situations 
where a single global transformation with only a few parameters can be used to relate the position of one image in a 
sequence with respect to neiqhbourinq images in a sequence. This will apply if the imaqe to be captured is planar (for 



